Static charged perfect fluid with Weyl-Majumdar relation 
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Static charged perfect fluid distributions have been studied. It is shown that if the norm of the 
timelike Killing vector and the electrostatic potential have the Weyl-Majumdar relation, then the 
background spatial metric is the space of constant curvature, and the field equations reduces to a 
single non-linear partial differential equation. Furthermore, if the linear equation of state for the 
fluid is assumed, then this equation becomes a Helmholtz equation on the space of constant curva- 
ture. Some explicit solutions are given. 
PACS number: 04.20. Jb 



I. INTRODUCTION 

In Newtonian theory, there are static equilibrium con- 
figurations of charged particles in which the gravitational 
attractions and the electric repulsions among the parti- 
cles exactly balance. In general relativity, it is well known 
that the Einstein equation has such a static solution cor- 
responding to these Newtonian situations, namely the 
Majumdar-Papapetrou solution P,pj describing an equi- 
librium state of many charged black holes. Q 

Majumdar studied the static Einstein-Maxwell sys- 
tem that is described by the static metric 



g^ fdt®dt- f^^h 



and the vector potential 



A = -^pdt, 



(1) 



(2) 



subject to the coupled Einstein and source-free Maxwell 
equations. The function / < is the square of the time- 
like Killing vector field dt, h is the Riemannian 3- metric, 
and the function tp denotes the electrostatic potential. 
Majumdar B proved that if there is a functional rela- 
tionship between / and ip, this must be of the form 



/ = -v^' 



(3) 



with a constant c, generalizing Weyl's result R] in the 
case of axial symmetry. Moreover, he showed that if 



./ = -V'' 



(4) 



then the 3-metric /i is a flat metric, i.e., h — dx^ + dy^ + 
dz^, and the field equations reduce to a single Laplace 
equation in the 3-dimensional flat space. Solutions of 
this type were studied by Papapetrou g] also. Later, 
Hartle and Hawking [Q investigated the global structure 
of the Majumdar-Papapetrou solution, and they showed 
that this can be interpreted as describing an equilibrium 
state of extremely highly charged black holes. 

For a system of charged dust in Newtonian theory, a 
similar equilibrium between gravitational and Coulomb 
forces is possible, whenever the mass density fi and charge 
density a of the charged dust satisfy 



(T ~ ±/x. 



(5) 



Das g studied the corresponding situation in general 
relativity, and he showed that the condition d|), as (H), 
implies the Weyl-Majumdar relation (^, and in this case 
the fleld equations reduce to a single non-linear equa- 
tion; i.e., the balance condition (0) makes sense in pre- 
cisely the same form in general relativity. Further discus- 
sion concerning the generality of the assumption (H) can 
be found in De and Raychaudhuri. H Explicit solutions 
for this s_^tem have been found by many authors, e.g., 
Refs. ||9|-|l^ with spherical symmetry, and a recent paper 
p3[ without spatial symmetry. 

The balance equation and the Weyl-Majumdar rela- 
tionship are closely related, and moreover, it seems that 
the force balance condition greatly simplifies the Einstein 
equation. The most typical examples are the Weyl and 
Majumdar-Papapetrou metrics mentioned above, and 
the Giirses metric with static charged-dust source. |l^] 
In Newtonian theory, the gravitational and electrostatic 
forces acting on an element of charged dust must be par- 
allel. This is also the case in Einstein theory; namely, 
the Bianchi identity (Euler equation) implies that there 
is a functional relationship between the gravitational and 
electrostatic potentials |§] (not necessarily in the form of 
(^ or (Q)). Such an analogy between the Newtonian 
and Einstein theory might hold even in the presence of a 
third force, e.g. a magnetic field or a pressure gradient. 
A configuration of a magnetic field parallel to the electric 
field is easily obtained by a duality rotation of pure elec- 
tric field. However, there will be a magnetic monopole in 
this case. The presence of non-parallel magnetic and elec- 
tric fields generates a Poynting fiux, contributing to the 
energy-flux component of the stress-energy tensor, which 
is a characteristic of a stationary field rather than static 
one. On the other hand, the (isotropic) matter pressure 
itself does not generate any energy (heat) fluxes. Rather, 
it effectively contributes only to the energy-density and 
momentum-flux parts of the stress-energy tensor by def- 
inition, and one may expect that there is a large class 
of static conflgurations generated by a charged perfect- 
fluid. Such a configuration may be considered as the 
interior metric of multiple Majumdar-Papapetrou black- 
holes. Since the Majumdar-Papapetrou metrics describe 



multiple (extremely highly) charged black- holes, they in- 
evitably have timelike singularities in each horizon, so 
that consideration of the regular interior metric might 
be interesting. Gautreau and Hoffman ||^ generalized 
the balance equation in the presence of matter pressure, 
assuming Eq. (j^), and showed 



IBttj 



* d* d^, 



a = ±{n + 3p) 1 



?/;2 



(6) 



where p is the isotropic pressure of the fluid, though they 
did not consider the rest of the field equations. The ex- 
istence of such a simple balance equation seems to im- 
ply that there are many static configurations. In fact, 
Guilfoyle [|l^ recently studied this system of a charged 
perfect fluid subject to (||) and obtained a large class of 
algebraically special (type D) metrics. 

In this paper, we investigate a system of a 
pure-electrically charged perfect fluid with the Weyl- 
Majumdar relation from a general point of view, and 
we show that there exists a large class of exact solu- 
tions, which can be algebraically general (type I) and 
need not have any spatial symmetries. As shown be- 
low, if the Weyl-Majumdar relation (Q) is assumed, then 
the 3-metric h must be that of a space of constant curva- 
ture, and the field equation reduces to a single non-linear 
equation. Moreover, if the fluid obeys a linear equation 
of state {p — const, x fi), then the field equation becomes 
a Helmholtz equation on the space of constant curva- 
ture (j|ll|). Therefore, solutions are expressed in terms 
of eigenfunctions of th e La place operator. Some explicit 
solutions are given in §111, 



II. FIELD EQUATIONS 

The metric of static spacetimes can be expressed in the 
form 



1 

g = -—rdt®dt + (l)'^h, 



(7) 



where neither the function cj) > nor the Riemannian 
3-metric h depends on the coordinate i; i.e., the vector 
field £, = dt satisfies the Killing equation 



A.9 - 0. 



(8) 



The Killing vector field <^ of a static spacetime is hyper- 
surface orthogonal, or twist-free, and thus it satisfies the 
Frobenius condition 



C A dC = 0, 



(9) 



where ^ — —(j)^'^dt is the Killing 1-form. (Here and in 
what follows, we use the same symbol for vector fields 
and their corresponding 1-forms.) 

The time-space components of the Ricci tensor are ob- 
tained from the equation (see e.g., Ref. [nsl) 



where 



Sttj = -ffic(^,(9^)da;^. 



(10) 



fll) 



and * represents the Hodge operator. The Lichnerowicz 
theorem states that j is proportional to ^ in the static 
case. [Q This can be shown by 

Q = d*{ihd£,) = dL^* di 
= [C^ — L^d) * d^ — 16ttl^ * j 
= -16Tr^Aj, (12) 

where ij denotes the interior product with respect to ^, 
and the identity £5 = di^ + i^d has been used. 

The time-time component of the Ricci tensor is ob- 
tained by calculating the Laplacian of / = — (/)~^: 

d* df — d* di,^^ = d * (£j — i^d)^ 

= d(^ A *dO = *(df ■.dO~^Ad*d^ 

= * 7 hl67r*5 '(C,j)- (13) 

Here the colon denotes the inner product of differential 
forms (i.e., a : (3 = (l/p!)a^^^2...^p/3'^^''^'"^p for p-forms 
a and (3). Therefore, we have 

Sttp = Ricit = 4>'^*d^d(j) + (j)-'^g-\d(t), d(j)). (14) 

The Ricci tensor of the 3-metric h encodes the remain- 
ing components of the Ricci tensor, and with a direct 
calculation, we obtain the expression pOl 

Riclh] = Ric ~ 8TT(t)^(pg -f 2^ ® j) + 2 "^"^^ . (15) 

Next, consider the static electro-magnetic field F = 
dA, C^F — 0. Here we assume that only the electric field 
exists for static observers, which implies 



CAF = 0. 



(16) 



In the source-free Einstein-Maxwell static case, this as- 
sumption is justified (up to duality rotation) by the Ein- 
stein equation, while in the charged perfect fluid case, 
this condition is equivalent to imposing the condition 
that the fluid element as well as the gravitational field 
is static. Moreover, there is no freedom of the duality ro- 
tation in this case. Here, the Maxwell equation dF — 
becomes 



di^F = {C^ - L^d)F = 0, 



(17) 



while the other Maxwell equation, d*F = 4TTa*u, reduces 
to 

d*L^F = -d{£, A *F) == - * (dC : F) + Ana * g^^{^, u) 
g-\L^F,d4>) 



= -2* 



+ A^a*g-\i,u), (18) 



where u and a denote the velocity and the charge den- 
sity of the fluid, respectively. Equation (|l^) implies that 
locally there exists a real function (electrostatic poten- 
tial) -0 such that dip = i^F holds. Then the full Maxwell 
equations reduce to a single Poisson equation, 

^d*di^^2?^^^^^P^-AT:ag-\S„u). (19) 



On the other hand, the stress-energy tensor of the 
Maxwell field Tp is 

SttT^ (a^, 9,) = g-\i9^F, is^F) + g-\ia^ * F, ta„ * F) 

(20) 

In pure electric cases, jp — —TF{^,d^)dx^ is propor- 
tional to ^. [Q Thus, the energy-flux components of the 
stress-energy tensor vanish, 

^AjF=0, (21) 

while the energy-density component becomes 

STrpf = SttT^ (C, a = 9'\diP, diP). (22) 

In terms of jf and pp, Eq. (M) can be written as 

SttTf = 8Tr(f)^{pFg + 2^® Jf) - 20^ dV' «i dip. (23) 

Next, consider a perfect fluid. Its stress-energy tensor 
has the form 



Tm = (m + P)u ®u+pg, 
or equivalently 



(24) 



Tm ^T„i - -tr(T„)g = {p + p)u (g) u + -{p~p)g. (25) 

The fluid velocity u must be proportional to ^; i.e., u = 
(/)^. (This is a consequence of the time-space components 
of the Einstein equation.) In analogy to the electric field, 
we define 



1 



jm = -T^m(C, d^)dx^ - -ifl + 3p)C 



and 



Pm =2i„(C>0 



(m + 3p) 
2(/)2 ■ 



The equation of hydrostatic equilibrium (Euler equation) 
is obtained from the Bianchi identity for the stress-energy 
tensor V • (Tp + T,n) = 0, which gives 



dln( 



1 



-{dp + acpdip). 



(28) 



p + p 
Now we consider the Einstein equation 

RtC^87T{TF+f,n). (29) 



The time-space component equations of the field equation 
^ /\ {j ~ 3f — jm) = are identically satisfied, while the 
time-time component p ^ pp + Pm leads to 



* d * d4) = — 



g-\dcp,d4>) 



(p-'g-\dip, dip) + Att(P{h + 3p). 
(30) 



From Eqs. (P_5f), ( p3| ) and (Eq), we obtain the remaining 
equation, 

Ric[h] + Wncp^ph = 2(p-^d(p ® dcp ~ icp^dip ® dip. (31) 



Our task is to solve Eqs. (|T|), (|28|), (|^) and (H). How- 
ever, we need additional conditions, such as an equation 
of state, to make this system deterministic. 

To solve Eq. (^, we assume the Weyl-Majumdar re- 
lation (Q), or 

iP + e(p-^=Q, (e = ±l). (32) 

Then the r.h.s. of Eq. ( pl| ) vanishes, and the 3- 
dimensional manifold with metric h becomes an Einstein 
space; i.e., the Ricci tensor has only a trace part (scalar 
curvature). Since in three or higher dimensional Einstein 
space, the scalar curvature must be constant (which can 
be seen from the contracted Bianchi identity), the matter 
pressure p is proportional to 0~^ with a constant ratio. 
Without loss of generality, we may take p = — k/Stt^^ 
(k = 0, ±1). Since a 3-dimensional Einstein space is a 
space of constant curvature (namely, the metric as well 
as the curvature is charactrized by a constant scalar cur- 
vature), we have 



dx^ + dy^ + dz'^ 
h = ; 

[l+iK(x2+2/2+2;2)]' 

Then, Eqs. {^ and (M) reduce to 



and 



(A-f A)(/) = 



a = €{p + 3p), 



(33) 



(34) 



(35) 



where A denotes the Laplacian with respect to h, and 
(26) A — A7T(p'^(u + 3p). It is easily checked that the Euler 



(27) 



equation (Eq) is consistent with (B5|). 

Equation (|34| ) is non-linear. However, it can be made 
linear if we assume a linear (isothermal-type) equation 
of state for the fiuid, i.e. p = const x p(= const x 0^^). 
Then, A becomes constant, and Eq. ( p4| ) is simply a 
Helmholtz equation in the space of constant curvature, 
whose solutions are known (see, e.g., Ref. [pO[). 



III. EXAMPLES 

Here we study some explicit solutions of Eq. (pj) . We 
mainly consider the spherically symmetric solutions as 
simple examples. We assume the dominant energy con- 
dition, which requires p > \p\. 



A. Positive pressure (k — —1) 

In this case, the dominant energy condition requires 
A > 2. In the spherical polar coordinate system {r, i?, (p}, 
the 3-metric h can be written as 

h = r]{rf[dr^ + r^{dd^ + sin^ t?V)], (36) 

where 77 = (1 - r^/4)"^ (0 < r < 2). Equation (|4|) for 
the function (/)(r) becomes 



1 d 9 (i , 

-r]r^ — (j) + X(j) = 0. 



The general solution of this equation is 

2 + 7 



,l-rV4 
1 = I sm 



(A -1)5 In 



S 



(37) 



(38) 



where / > and < 6 < tt are integration constants. 
In general, the set on which = is a singularity, since, 
for example, the field invariant F : F = —(j>~'^g'~-^{d(p, dcj)) 
blows up there. Hence this solution must have a singular- 
ity, and the situation is the same even without spherical 
symmetry. The locations of the singularities are given by 



2tanh 



(m: — S) 
2(A-1)^ 



(n = 0,±l,±2, 



(39) 



at which the area radius R — \4'rir\ vanishes, so that 
these singularities are actually point-like, and each con- 
nected component of the spacetime has in general the 
topology S^ minus two points. (Without spherical sym- 
metry, the singularity might be line-like or surface-like, 
and the spacetime topology would be S^ minus several 
such sets.) 

Next, let us consider the solution (|3|) in the neighbor- 
hood of r = 0. We may assume that < S < tt/2, since 
the transformation r 1— > —r corresponds to 6 *—> n ~ S. 
Then, let us consider the region r < 0. Whenever <5 7^ 0, 
this solution diverges at r = like IsinS/r. However, 
this is only a coordinate singularity. The hypersurface 
r = is a null hypersurface, on which the Killing vector 
becomes null, and it is generated by spheres of radius 
R — / sin cr. In fact, whenever 6 ^ 0, 7r/2, the region 
r < is extended analytically to the region r > across 
this null hypersurface. To see this, let us consider the 
coordinate transformation 



V ^t + 



(p rjdr. 



(40) 



Then the metric transforms to 


(p^ 1 — r^/4 


+ f sin^ 


1 2 + r' 



{dd^ +sin^Mip^), (41) 



which is manifestly non-singular at r = 0. 



Moreover, we can consider a class C^ extension of the 
region r < to the exterior field of a black-hole space- 
time. The Misner-Sharp mass mn (which coincides with 
the Hawking mass or Kodama mass in the case of spher- 
ical symmetry) with respect to the surface r = is 



/ sin (5 



rriH 



(42) 



On the other hand, the electric charge inside a two- 
surface S can be evaluated as 



In particular, when a surface defined by r 
chosen as S, we find 

r2 d\(j>\ 



(43) 



const is 



1 - r2/4 dr 



(44) 



where 'inside' corresponds to smaller r, so that the elec- 
tric charge inside the surface r = is 



qh = d sin S = 2eniH- 



(45) 



This relationship is identical to that satisfied by the 
extreme Reissner-Nordstrom spacetime on the horizon. 
Furthermore, the mass density and matter pressure of 
the fluid vanish on the surface r — Q. These evaluations 
suggest that the region r < with the metric charac- 
terized by the Eq. (JS^) can be matched at r = to the 
exterior field with the extreme Reissner-Nordstrom met- 
ric. By rescaling the coordinates according to 



l{\ - 1)2 cos 5 r, i = 



t 



1{X — 1)2 cos (5 



the metric transforms to 


g = "l^' -1- 02^2 [^~2 _^ ~2 (^^2 ^ gj^2 ^^^ 
02 


where 


1 ^ ■ 
(p = -TT sm 

rr] 


2l{X-l)icos6-r 



and 



V 



1 



4?2(A-l)cos2(5 

Note that behaves near f — like 
/sin (5 



= 1 



0{f 



Next, we introduce the advanced time coordinate 



t + 



Z2- J- 



r)df. 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



Then the metric ( p7| ) becomes 



9 = 



+ 2f]dvdf + ^'^yff^idiP + sin^ Odip^), (52) 



while the metric of the extreme Reissner-Nordstrom has 
the form 



9rn 



1 + ^ ) dv^ + 2dvdr 
+ [l + ^)%\dd^+B\n^dd^''). (53) 



If the mass parameter (ADM mass) of the metric (m) is 

TO = Zsin(^, (54) 

then these two metrics are matched on the hypersurface 
f = 0, and it is easily seen from ( |50| ) that this extension 
is class C^. In this case, the exterior metric ( p3[ ) and 
the interior metric (52) correspond to f > and f < 0, 
respectively. We can similarly verify that the region r > 
for ( p8[ ) can be matched to the m = I sin 5 interior metric 
(i.e., the metric (^) for to, — —/sin (5 < 0, f > 0), and if 
I sin (5 < 0, then the regions r < and r > for (BS) can 
be matched to the interior (m = — |/sin(5|) and exterior 
(to = |Zsin5|) metric (|53|), respectively. 

Of course, we may also obtain solutions without spher- 
ical symmetry using the method of multipole expan- 
sion. Instead, here we give an example by superposing 
monopole solutions with different centers {xi,yi, Zi): 



E'.^ 



r2/4 



sm 



(A-l)Mn: 



where 



Tj = [{x - x,f + {y~ y^)^ + {y - Vif 



, (55) 



(56) 



Here, the quantities li and 5i are constants. Each hole 
Ti = Q can be matched to the Reissner-Nordstrom metric 
with m = k sin 6i . 



B. Dust K = 

For completeness, we briefly discuss this special case. 
We only give some explicit solutions here, since the prop- 
erties of these solutions are quite similar to those in the 
positive pressure case. Recently, Giirses |1^ also found 
this class. The field equation is the ordinary Helmholtz 
equation, 



dx'^ 









+ AU = 0, 



and the energy condition requires A > 0. 
spherically symmetric solution is 



(57) 



The general 



sin(A2r + S) 



(58) 



where / > and < 6 < tt are integration constants. 
As in the positive pressure case, the solution vibrates 
around (/> = in general, so that there always exists 
a singularity. The solution is matched to the exterior 
Reissner-Nordstrom metric with to = I sin 5 at r = 
with class C^ differentiability. Moreover, in the case of 
spherically symmetric static dust solutions, the metric 
can be matched to an extreme Reissner-Nordstrom met- 
ric on any symmetric sphere by virtue of the absence of 
matter pressure. 

The solution without spatial symmetry can be ex- 
pressed in terms of the spherical harmonics and the 
spherical Bessel functions. We simply give a multi-center 
solution 



E^ 



sin(A2ri -I- S) 



(59) 



For more information, see Rcf. llJ 



C. Negative pressure (k = 1) 

In this case the dominant energy condition requires 
A > —1. Consider the solution of Eq. (|3^) with 77 = 
(1 + r^/4)~^. The general solution is 



.1 



l- 



-r74 
sin 



(1 -I- A)2 arcsin 



H-r2/4 



+ S 



(60) 



where I > and < S < n are integration constants, and 
the arcsin above should be regarded as 



Arcsin 



arcsin 



H-r2/4 



Arcsin 



l+r^/4 



Arcsin 



l+r2/4 



(r < -2) 
(-2 < r < 2) 

(r>2) 

(61) 



in terms of the Arcsin (which takes the principal value). 
Unlike the case of positive pressure or dust, these solu- 
tions need not have any singularity in the appropriate 
range of A. In fact, if -1 < A < -3/4 {\p\ < ^ < 3/2|p|) 
and (1 -I- A)2 7r < S < tt — (1 + A)2 7r, then (f> gives a 
regular solution for all r, and ii S — and — 1 < A < 
{\p\ < M 5: 3|p|), then it gives a regular solution for r > 0. 
As in the other cases, we can match the solution to 
the exterior metric (f > 0) of the extreme Reissner- 
Nordstrom spacetime (^3|) at r = whenever S ^ 
0, 7r/2, TT. Note that we may assume < i5 < 7r/2, since S 
is transformed to tt — (5 by r 1— > — r. Let us consider the 
coordinate transformation 



f = /(H- A)2 cos 5 r, < = 



t 



l{l + X)2 cosS 



(62) 



Then the metric can be written 






^f2 
^2 


, (63) 


where 




(p = — sm 
ffj 


(] 


1 TT] 


, (64) 


n^ A. J ^ cLi Colli ^ n^ u 

1(1 + \)^ cosS 


and 




^ 


[l 1 '' 1 


-1 


(65) 


V ~ 


4/2(i + A)cos2,5 


Furthermore, 


the transformation 







V = t + 



^Tjdr 



(66) 



gives a metric in the form 

g = -}-^dv^ + 2fidvdf + 4>^fi^r^{d^^ + sin^ M(p^). (67) 

The region f < with this metric can be matched to the 
region r > with the metric (m) with 

m — I sin (5. (68) 

In this case also, the extension is class C^, since 

/ sin (5 



1 



+ 0(P), r7=l + 0(P). 



(69) 



Unlike the case of a positive pressure fluid, the metric 
( pOl ) or (|6^) can be regular for all r < (f < 0) if (5 > 
(1 + A)2 7r is satisfied. In this case, the surface r — — oo 
is neither infinity nor a singularity. In fact, it can be 
shown that this surface is a null hypersurface and that 
the spacetime can be extend beyond r = — oo, unless 
(5 = (1 + A)2 7r, in which case r — — oo simply corresponds 
to a regular center. To see this, consider the inverse of 
the radial coordinate, 

r = 4/r. (70) 

In terms of this coordinate, the metric transforms to 

de 



9 



where 



+ (jj^f [df2 + f2(dt?2 _^ sin2 ^^^2)j ^ (-7-^-) 



l+f2/4 

/ sin 



(1 + A)2 arcsin 



l + f2/4 



+7r + (1 + A)27r-(5 



and 



1 + -P 
4 



(72) 



(73) 



Hence, the effect of this transformation is simply 

S^7r+{l + X)K-6, (74) 

so that the metric can be analytically extended across 
the null hypersurface f = in a manner similar explained 
above, and it can also be matched to the interior region 
(f < 0) with the Reissner-Nordstrom metric (pSh with 
TO = /sin[(5 — (1 + A) 2 7r]. 

Finally, we obtain a multi-center solution in an obvious 
way; 



E^- 



r2/4 



sm 



(1 + A)2 arcsin 



l + r2/4 



+ S, 



(75) 



IV. CONCLUSIONS 

We have investigated the static field of a charged per- 
fect fluid, and shown that only two assumptions, the 
Weyl-Majumdar relation (m and a linear equation of 
state, allow us to reduce the field equations to a sin- 
gle linear equation. Thus, any solution of this type can 
be described by an eigenfunction of the Laplacian in the 
3-dimensional space of constant curvature, and the eigen- 
value determines the equation of state, i.e., the ratio of 
the energy density to pressure. These solutions have 
some desirable properties. For example, they obey a 
simple equation of state, as required, neither the mass 
density nor the matter pressure changes sign, and they 
need not have any spatial symmetry. 

The mechanism responsible for this simplification of 
field equations may be understood in terms of force bal- 
ance among the gravitational force, the Coulomb force 
and the pressure gradient. The function cj)^"^ is the norm 
of the static KiUing vector, and thus it can be interpreted 
as the gravitational potential. Then Eq. (||) gives a func- 
tional relationship between the gravitational and electro- 
static potential, which implies that the gravitational and 
Coulomb forces are parallel. This property characterizes 
the Weyl ||] and Majumdar-Papapetrou [||j|] solutions. 
One main point of this paper is that the pressure gradient 
automatically becomes parallel to both the gravitational 
and Coulomb forces as seen from the Einstein equation, 
even in the presence of a perfect fluid, which is implicit in 
the spherically symmetric case. Thus, the balance equa- 
tion among the three relevant forces, which is originally a 
vector equation, becomes a scalar equation (p5|), explain- 
ing why a large class of inhomogeneous metrics are ob- 
tained. Though we have discussed only the situatioin in 
which equipotential surfaces of the gravitational poten- 
tial, the electrostatic potential, and the matter pressure 
coincide, it is not a necessary condition for the static con- 
figuration; there would be a situation in which the vector 
balance equation holds at each point. However, such a 



spacetime might have some spatial symmetry, such as 
axial or cylindrical symmetry, or the equation of state 
might depend on spacetime points. The static configura- 
tion considered here may be regarded as a final state of a 
non-rotating object composed of charged fiuid. Through 
the realistic contraction of such an object with a given 
equation of state, the scalar balance equation might be 
satisfied in the final equilibrium state, even if it is not 
spherically symmetric. 

Another characterization of the solution can be given 
in terms of the Ernst potentials pll. In analogy to 
the stationary Einstein-Mawell system, we can define the 
Ernst potential by £: = l/(jj^ - V'^ Then Eq. @ can be 
replaced by the Poisson equations for £. However, the 
assumption (y) implies that £ = 0, so that this equa- 
tion reduces to an algebraic equation, which is just the 
balance equation (p5|). 

We have also shown that solutions corresponding to 
Eqs. (|38|), (|5q ) and ( pOf ) can be matched to the extreme 
Reissner-Nordstrom metric at the horizon. However, this 
interior solution has a central singularity unless the mat- 
ter pressure is negative. In the negative pressure case, the 
interior solution d7q) can be regular for an appropriate 
range of parameters, and then r ^ — oo does not rep- 
resent conformal infinity but, rather, represents another 
'hole' connected by another interior metric, or regular 
center. Since the matching conditions are determined 
by the behavior near the horizon, solutions ([75[), (|59| ) 
and ( p5| ) could be matched to an appropriate Majumdar- 
Papapetrou metric ^-^1 at r^ = in a manner similar to 
that described here. 

The existence of a regular interior metric has some con- 
nection with the Buchdahl theorem. |23] Roughly speak- 
ing, the Buchdahl thorem states that a relativistic star 
composed of perfect fluid with mass M and radius R can 
exist only if i? > (9/4)M. The charged generalization of 
the Buchdahl theorem has been proved by Yu and Liu. 
p4| They showed that the radius of a star of charged fluid 
has a similar lower bound if the total charge of the star 
does not exceed the extreme limit Q — M . In the ex- 
treme case, there is a series of interior solutions in which 
the surface of the star can be arbitrarily close to a black- 
hole horizon, |£^,^ showing that inf i? = M. On the 
other hand, our example shows that there actually exists 
an i? = M star. 
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